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Abstract 

We suggest that the lithiation of pristine SnO forms a layered LixO structure while the 
expelled tin atoms agglomerate into ’surface’ planes separating the LixO layers. The proposed 
lithiation model widely differs from the common assumption that tin segregates into nano¬ 
clusters embedded in the lithia matrix. With this model we are able to account for the various 
tin bonds that are seen experimentally and explain the three volume expansion phases that 
occur when SnO undergoes lithiation: (i) at low concentrations Li behaves as an intercalated 
species inducing small volume increases; (ii) for intermediate concentrations SnO transforms 
into lithia causing a large expansion; (hi) finally, as the Li concentration further increases a 
saturation of the lithia takes place until a layered Li20 is formed. A moderate volume expan¬ 
sion results from this last process. We also report a ’zipper’ nucleation mechanism that could 
provide the seed for the transformation from tin oxide to lithium oxide. 

KEYWORDS: Lithium ion batteries, Anode, Tin oxide, Lithiation, Nueleation, Density fune- 
tional theory 

1 Introduction 

Tin based compounds represent a promising category of materials to be used as the anode of 
lithium-ion batteries (LIBs).®With an observed capacity of either 4.4 (Li22Sn5)l21or 4.25 (Lii 7 Sn 4 )l^ 
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Li per Sn atom, they outperform the maximal theoretieal energy density of graphite eleetrodes by 
a faetor greater than two.® The interest in Sn-based systems was triggered by the work of Idota 
et a/.® who applied an amorphous tin eomposite oxide as the anode material. Several other re- 
seareh groups have later investigated similar eompounds in the eontext of LIBs.®® However, a 
eommereially viable solution remains to be determined sinee Sn-based systems tend to suffer from 
an unaeeeptably fast fade in eapaeity. This happens in spite of the various attempts to optimize the 
ehemieal eompositionl^HIIl and struetural features^^^^ of the anode, showing that a better under¬ 
standing of the atomie strueture of the aetive materials is required for improved designs. A reeent 
review by Park et a/.® lists the most relevant advanees made in Sn-based eleetrodes and their main 
eauses of degradation. The large volumetrie ehanges of Sn, SnO, and Sn02 during lithiation and 
delithiation is one of the main deterioration meehanism and among the hardest teehnologieal ehal- 
lenges to overeome. An example of the expansion taking plaee when SnO is lithiated is displayed 
in Figurewhieh reports the volumetrie ehange of LixOSn with respeet to the initial SnO erystal 
as a funetion of the Li load level.® At moderate Li eoneentrations, i.e. up to a ratio of two Li atoms 
per Sn atom, three distinet expansion regimes ean be identified from this experimental data, as 
diseussed later. 

In one of the first studies addressing the struetural ehanges oeeurring during the lithiation of 
Sn-based eompounds (SnO, Sn02, Li 2 Sn 03 , and SnSiOs) Courtney and Dahn® suggested that the 
initial tin oxide transforms into an amorphous lithia matrix, from whieh the Sn atoms are expelled 
before segregating into nano-elusters embedded within the forming matrix. This lithiation model 
relies on two phenomena. First, an irreversible proeess takes plaee during the first lithiation eyele, 
whieh ean be attributed to the formation of a stable lithia matrix. Seeondly, XRD measurements 
of the lithiated samples eontain signatures of Sn-Sn and Sn-Li bonds, potentially validating the 
proposed elustering of Sn atoms. In a later work the transformation towards a lithia was eonfirmed 
by Chouvin et al^ who supported this with results from Mossbauer speetroseopy. They revealed 
that the Sn atoms are redueed as the Li eoneentration level inereases. However, their data also 
showed that the forming lithia does not eorrespond to its bulk eounterpart sinee markers for various 
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unknown bonds were present. A similar eonelusion was reached by Sandu et a/.^^who also applied 
Mossbauer spectroscopy. Their results agian showed signatures of unidentified Sn bonds after the 
transformation into lithia had taken place; they named them ’exotic’ bonds. Furthermore they 
observed that a spectra consistent with SnO was present at low Li concentrations and therefore 
concluded that Li behaves as an intercalated species during the initial part of the lithiation process. 
As the Li load level increases they also found that parts of the SnO matrix could still accommodate 
Li atoms as intercalated species, while other parts started transforming into lithia. 

Another peculiarity of Sn-based systems was later reported by Zhong et a/.^and further inves¬ 
tigated by Nie et a/.^In their experiments they saw that the transformation into lithia occurs along 
stripes when a nanowire of Sn02 gets floated with lithium. Such a nucleation is very different from 
the usually assumed core-shell mechanism and remains to be understood. Although the crystalline 
structure of SnO and Sn02 are significantly different it is worth investigating whether Li preferably 
arranges itself along stripes in other Sn-based oxides. From the experimental results summarized 
above it appears that (i) three expansion regimes can be seen when SnO is lithiated,® (ii) Li atoms 
behave as intercalated species at low concentrations, (iii) ’exotic’ Sn bonds manifest themselves 
after the first lithiation cycle, and (iv) an abnormal nucleation mechanism occurs in Sn02. None of 
the existing lithiation models simultaneously captures all these phenomena and offers a consistent 
explanation for them. Moreover, the common models assuming a complete segregation of Sn and 
Li20 underestimates the experimentally determined volume expansion by 7.5%, clearly indicating 
that important effects are not taken into account. 

To address these issues we report here results from simulations investigating the structural 
evolution during the lithiation of pristine SnO. Starting from a SnO crystal, atomistic simulations 
are conducted to determine the volume expansion as the oxide transforms into lithia up until the 
point where a saturated Li20^ is formed. The final structure is layered and different from its 
crystalline counterpart and is therefore marked with a superscript L. We also demonstrate that 
rather than segregating into nano-clusters the Sn atoms remain evenly distributed between the 
forming Li 20 ^ layers so that they can be considered as ’surface’ planes on the appearing lithia. 
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This new model ehallenges the eurrent understanding of SnO lithiation as being a process where a 
Li20 matrix forms and the expelled Sn atoms segregate into nano-clusters. The paper is organized 
as follows: we start by describing the computational approach. Simulation results are listed in 
Section which is divided into two parts, one dealing with tin oxide and the other with lithia. 
Section discusses the obtained results and their consistency with experimental data. Finally, 
conclusions are drawn in Section [51 



Capacity Q [mAhg^] 


Figure 1: Relative volume expansion of LixOSn with respect to a SnO sample during the first 
lithiation cycle. The reported curve only represents the initial part of the full lithiation/delithiation 
cycle that is shown in the inset. The black line refers to the experimental results of Reference |8] 
The solid black square, which is also taken from Reference |8l represents the volume expansion 
for the case where Li20 and Sn completely segregate. Three regimes can be identified from the 
measured curve, which are indicated by the colored lines (i) red, intercalation (ii) blue, oxide 
transformation and (iii) green, saturation of the lithia. The crosses are results from our simulations. 
They are colored according to the expansion regime they correspond to. The dashed lines connect 
computed data points belonging to the same regime. 


2 Computational Details 

Two types of supercells are utilized to determine the structural changes of SnO as it lithiates. In 
the first series of simulations addressing low Li concentrations in SnO and the formation of the 
Lii/ 20 Sn crystal, the supercell consists of two layers, whereas it only contains a single layer for 
the calculations of LiOSn and Li20^Sn. The definition of a layer is given in Figure Both 
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Figure 2: Pristine SnO (color online). Yellow colored atoms represent tin and red oxygen. Only 
bonds shorter than 2.3 A are shown. It is know that tin oxide is a layered and ordered structure.^SI 
Here three such layers are shown and it appears that the (100) and (010) planes are identical. 
Each layer consists of two flat Sn planes and a slightly wiggling central plane of O atoms. Within 
each layer the Sn and O atoms connect to four neighbors of the opposite species. The center 
most layer—enclosed in the dashed rectangular box—forms one of the two layers included in the 
simulated supercell. 


supercell types exhibit a total of 64 SnO units (64 Sn and 64 O atoms), while the number of Li 
atoms varies from 0 to 128. 

The atomic interactions are accounted for by applying the VASF^ density functional theory 
(DFT) tool. Core electrons are described by the projector augmented wave (PAW) method.l^ All 
simulations are performed within the Perdew-Burke-Emzerhof (PBE) parametrization of the gen¬ 
eralized gradient approximation (GGA)^ to the exchange-correlation functional. The KPOINT 
sampling is 2x2x2 and the plane-waves are truncated at 500 eV. The electronic degrees of freedom 
are considered as converged when the change in energy decreases below 1 ■ 10^^ eV per self- 
consistent field (SCE) iteration. A Nose-Hover thermostat at either lOOOK or 500K is applied in 
certain cases to enable transitions between nearby local minima. Einally, all structures are relaxed 
in terms of their ionic degrees of freedom (atomic nuclei) and the size and shape of the super¬ 
cell through the conjugate gradients method as implemented in VASP. In this case convergence is 
reached when the force decreases below 1 ■ 10“^ eV/A. 


5 



















3 Results 


3.1 Tin oxide 

At low concentrations the inserted Li atoms have a preferenee for linear struetures within the SnO 
matrix, i.e. they are favorably arranged as stripes. The lowest energy configurations result when 
the Li atoms reside in a single Sn-O-Sn layer and neighbor Li atoms belong to alternating Sn 
’surfaces’. At this concentration level the Li atoms behave as intercalated speeies and are likely 
subject to diffusion within the plane spanned by a Sn-O-Sn layer. In the following two seetions 
these findings are explained in greater details. 

3.1.1 Pure oxide 

A pristine tin (II) oxide crystal is generated from a tetragonal PbO lattice in which the Pb atoms 
are replaced by Sn and relaxed.l^ The resulting bulk phase of SnO has a distinet layered Sn-O- 
Sn structure in the Z-direction, which clearly appears in Figure that shows three such layers. 
Eaeh of them ean be divided into three planes: a eentral plane of four-fold coordinated O atoms 
in tetrahedral sites is sandwiched between two ’surfaee’ planes of four-fold eoordinated Sn atoms 
that reside at the apex of a square pyramid. The bond length between Sn and O is 2.3 A, whieh 
agrees reasonably well with the experimental determined value of 2.21 A^^and a distance of 2.3 A 
is used as an the upper value to determine whether two atoms are eonnected by a bond. All the 
atoms have four nearest-neighbors and the resulting supereell is a tetragonal box with identical 
(100) and (010) planes. The distanee between two O planes is 4.9 A and the footprint per SnO 
unit in the (001) plane is 7.4 A^ for a total volume of 36.7 A^ per unit. Within the SnO structure 
the distance between the upper and lowermost Sn ’surface’ planes of two adjacent layers is 2.6 A, 
which exceeds the height of the layer itself. 

In the limit of a single Li atom added to the SnO matrix, the site preference is within one of 
the Sn planes. The lowest energy eonfiguration is found when the Li atom occupies the same site 
as a Sn atom belonging to the opposite ’surfaee’ separated by the central O plane. In other words. 
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the O plane aets as a mirror between the newly introdueed Li atom and an existing Sn element. In 
this eonfiguration the binding energy for the atom is -250 meV and it is not thermodynamieally 
favorable eompared to the eohesive energy of a Li atom in a bulk environment. When loeated in 
this site the Li atom does not eause any signifieant disturbanees of the SnO matrix and the bond 
length between neighbor Sn and O atoms remains equal to 2.3 A. The distanee from Li to O is 

2.2 A. The separation between Li and the elosest Sn atoms in the same plane is 2.8 A, while the 
distanee to its mirror Sn atom residing in the opposite plane is 2.5 A. These faets eombined with the 
relative large spaeing between the SnO layers elearly show that at low eoneentrations the Li atoms 
ean be eonsidered as interealated speeies within the SnO matrix. Furthermore, at this eoneentration 
level, Li atoms are likely subjeet to 2D diffusion within the (001) plane. 

3.1.2 Lithium at low concentrations 

To investigate how lithiation affeets the SnO matrix at low eoneentrations, numerous atornie eon- 
figurations are generated where the load level inereases from two to eight Li atoms. In the seareh 
for the lowest energy strueture, eonfigurations are sampled with the additional Li atoms arranged 
to (i) reside within the same or neighbor SnO layers, (ii) be nearest-neighbors or separated, and 
(iii) ordered in different patterns. Struetures with four or more Li atoms are generated by embed¬ 
ding Li into a low energy eonfiguration with less Li atoms. This teehnique allows for a more rapid 
eonvergenee towards the lowest energy struetures. To relax the stress eaused by the inserted Li 
atoms, the heat treatment mentioned in Seetion|^is applied. The determined lowest energy strue¬ 
tures for three of the investigated load levels are shown in Figure A eomplete list of the tested 
eonfigurations ean be found in the supporting information. 

When the seeond Li atom is inserted into the SnO matrix the ratio between Li and Sn beeomes 
2:64. The two Li atoms then have the possibility to arrange themselves in a dimer eonfiguration. 
A weak binding of 42 meV/atom is aehieved when they form a symmetrie strueture and reside 
in opposite Sn planes. By eomparing to the strueture with a single Li atom in the SnO matrix it 
appears that 292 meV is gained per atom and an energy barrier of 584 meV has to be surpassed for 
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the dimer to break. In this eonfiguration the distanee to the neighbor O atoms is redueed to 1.9 A 
instead of 2.2 A as for the single Li atom. The separation between the Li atoms is 2.6 A, 0.1 A 
shorter than the distanee between two O atoms within the same plane. Finally, the bond length to 
the elosest Sn atoms within the same plane slightly inereases from 2.8 to 2.9 A, while the distanee 
to the ones in the opposite plane within the same layer ehanges from 2.5 to 2.7 A. When the dimer 
forms the surrounding SnO environment is somewhat disturbed and a few bonds eonneeting Sn 
and O atoms are stretehed beyond the defined bond value, as ean be seen in Figure]^ The volume 
inerease is insignifieant as eompared to the pristine SnO. 

If an additional pair of atoms is inserted, a total of four Li atoms is embedded in the SnO 
matrix and a 4:64 ratio is reaehed. At this load level the lowest energy eonfiguration is a eluster 
with two atoms in eaeh of the Sn ’surfaees’ as shown in Figure [^(b). The eorresponding volume 
expansion is 1.03%. The binding energy now inereases to 266 meV/atom, whieh is the highest 
value per Li atom in any of the inspeeted low eoneentration struetures. The separation between 
the eentral pair of Li atoms is 2.7 A, whereas the distanee between the rightmost pair is 2.6 A, and 
2.4 A for the leftmost one. The elosest Sn atoms are situated in the opposite plane at a distanee 
ranging from 2.7 to 2.9 A. The shortest distanee from Li to a neighbor O atom is in all oases olose 
to 1.9 A. It is noteworthy that the Li-0 bond length has a very narrow spread. In faet, this value 
eorresponds to the Li-0 bond length in the fully saturated Li20^ layer, whieh is deseribed in a 
later seotion and stands out as a strongly desired struotural feature. The four-atom eluster laoks 
symmetry beoause it seems impossible to oonform to the surrounding SnO matrix (bond length of 
2.3 A), while preserving the bonds to the O atoms. 

When six Li atoms are present in the oxide matrix a ratio of 6:64 is reaehed. The relaxation of 
the atomie strueture now leads to a replaeement of one Sn atom in the oxide matrix by one of the 
inserted Li atoms. As a eonsequenee the Sn atom is pushed outwards and away from the eentral 
plane when the Li atom situated in opposite plane enters the O plane. This transition spontaneously 
oeeurs indieating that no energy barrier needs to be surpassed for lithia to nueleate. After the 
transformation the binding energy per Li atom slightly reduees as eompared to the eonfiguration 
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containing a cluster of four Li atoms and is equal to 240 meV/atom. The volume expansion reaehes 
1.55%. In the relaxed strueture the distanees between neighbor Li atoms are 2.5 A in three of the 
five eases, whereas for the last two the distanee inereases to 2.9 and 3.1 A. The latter length 
eorresponds to the ease where a Li atom is within the eentral O plane. Again the shortest bond 
from Li to a neighbor O atom is 1.9 A and the elosets Sn atoms are within a range of 2.7 to 3.0 A. 

Under the assumption that linear arrangements of Li atoms remain at higher eoneentrations, 
the observed nueleation meehanism will aet as a ’zipper’ that joins the Li atoms and form a stripe 
of lithium oxide. The arrows in Figure [^(e) indieate the direetions of propagation for the Li atoms 
as they enter into the O plane. It is also important to note that during this transformation the Sn 
atoms will evenly be expelled above and below the emerging stripe of lithia. 

2 Li atoms, 42 meV/Li 4 Li atoms, 266 meV/Li 6 Li atoms, 240 meV/Li 



Figure 3: Interealation of Li in pristine SnO (eolor online). The same eolor eode as in Figure is 
used with the addition of blue for lithium, a) A Li dimer with an atom in the Sn planes above and 
below the eentral O plane is the lowest energy eonfiguration for two Li atoms, b) Four Li atoms 
form a eluster with two atoms in the planes above and below the eentral O plane. This results 
in the strongest binding energy per Li atom among all the samples where Li is eonsidered as an 
interealated speeies. e) Nueleation and ’zipper’ meehanism for the transformation from tin oxide 
to lithia. The arrows mark the direetion along whieh Li atoms should move to form a stripe of 
lithia. It is important to note that the Sn atoms evenly move to the ’surfaee’ sites situated above 
and below the forming LiO. 


3.2 Lithia 

At moderate Li eoneentrations a sudden aeeeleration of the volume expansion oeeurs, whieh is 
explained by the transition from SnO to LiO. We will not go into the details of this proeess, but 
rather foeus on the resulting Lii/20Sn, LiOSn, and Li20^ eonfigurations. Sinee these struetures 
are fully transformed oxides without any traees of interealated Li, the eomputed volumes are upper 
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Figure 4: Relaxed Lii/20Sn sample (color online). The same color code as used before. The 
dashed box encloses the simulated supercell. The blue lines are imposed to underline the striped 
structure of the Li atoms within the transformed oxide. As Li enters sites within the central O 
plane Sn atoms are pushed out. This causes a significant change in the bond configuration. For 
all the Sn atoms situated above the central O plane the coordination is reduced to three. Half of 
the Sn atoms below the central plane are now only bonded to two neighbor O atoms. The layered 
structure of the initial tin oxide remains, however, the supercell is skewed in the (001) plane and 
along the Z-direction. 


bounds for the corresponding capacities, which is also apparent in Figure[^ Despite the misaligned 
’volume expansion vs capacity’ these configurations are investigated to determine whether the 
layered structure of the initial SnO persists and if the expelled Sn atoms are still evenly distributed 
above and below the central plane. 


3.2.1 Lii/20Sn 

To gain insights into the structural changes taking place at moderate Li concentrations an atomic 
configuration with a 1:2 ratio between Li and Sn is constructed. The structure is made from a 
pristine SnO crystal where lines of Li atoms are inserted into the Sn planes above and below the 


central O plane. Based on the findings at low Li concentrations in Section 3.1.2 the Li atoms are 


placed at mirrored Sn sites. Every second row is left empty. Of all the simulated samples Lii/20Sn 
is the one that imposed the greatest challenge with several calculations resulting in layered, but 
disordered structures. However, it turns out that the structure with the lowest energy exhibits a 
layered and ordered arrangement, as shown in Figure It should be noted that the Li atoms 
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form linear chains, which resemble the stripes experimentally seen during the lithiation of Sn 02 .^ 
Although the materials are different (SnO vs. Sn02) linear chains of Li also appear in SnO. This 
could be a typical feature of Sn-based oxides. 

Compared to the initial SnO crystal, a drastic evolution is seen when Li atoms enter sites situ¬ 
ated within the central O plane, which causes a transformation from tin oxide to lithia. Following 
the modifications of the central plane, Sn atoms are pushed outwards and the majority of them 
share only three bonds with the neighbor O atoms. This should be compared to the four-fold co¬ 
ordinated sites within the pristine SnO crystal. The remaining Sn atoms located below the central 
plane are pushed even further outwards to sites where only a single bond to O remains. The in¬ 
serted Li atoms occupy sites where they are connected to three nearest-neighbor O atoms. All the 
O atoms keep their four-fold coordination and are located in tetrahedral sites with either two bonds 
to Li (1.9 A) and Sn (2.1 A), or one bond to Li (2.0 A) and three to Sn (2.2 A). 

The expulsion of Sn atoms and the oxide transformation increases the volume of the structure 
by 23.14%, as reported in Figure It also modifies the symmetry of the supercell from orthog¬ 
onal to triclinic. Despite these significant changes it is important to realize that the layered and 
ordered structure persists. Furthermore, the segregated Sn atoms do not seem to assemble into 
nano-clusters, but rather to comply with the planar structure of the SnO matrix. 

3.2.2 LiOSn 

As a next step LiOSn is generated by decorating all sites, within the Sn planes of the pristine SnO 
crystal with Li atoms. During the structural relaxation and energy minimization the intrusion of 
Li atoms continues and all the O atoms in the central plane are replaced by Li. The result of this 
transformation is depicted in Figure that represents a family of low energy configurations with 
some recombination in the interfacing Sn ’surface’ layers. Another example is reported in the 
supplementary material. In the selected configuration the central plane is occupied by four-fold 
coordinated Li atoms that reside at tetrahedral sites with bonds to O of length 1.9 A. The O atoms 
are now located between the Li and Sn atoms in either a five-fold coordinated square pyramid with 
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four bonds to Li in the base and a single bond to Sn in the apex (2.0 A) or a six-fold eoordinated 
site in a truneated eube with four Li atoms in the the base and two Sn bonds of length 2.2 A in the 
top diagonal. The Sn atoms are expelled from the eentral plane and are now residing as ’surfaee’ 
planes on the emerging LiO slab. By earefully inspeeting this strueture it turns out that the Sn 
plane above and below the lithia are slightly different. The Sn atoms on the uppermost ’surfaee’ 
reside at Li sites, whereas O sites are oeeupied in the plane below. Sn atoms originally belonging 
to two different neighbor layers are now loeated direetly above eaeh other in the (100) and (010) 
planes but the Sn atoms in the two layers are shifted in the (001) plane. Due to the ehosen eut off 
when plotting bonds the Sn atoms at Li sites (above) appear to have two bonds, whereas the ones 
at O sites (below) only exhibit a single bond. 

As eompared to the SnO referenee erystal the volume expansion of LiOSn is 29.42%, while it 
was 23.14% for Lii/20Sn. This moderate ehange indieates that the lithiation of SnO has reaehed 
a stable state, while the phase ehange (SnO —LiO) and a huge volume inerease oeeurred during 


the previous step deseribed in Seetion 3.2.1 Again, the ordered and layered strueture of the initial 
SnO sample remains. The sides of the supereell are skewed, as for the Lii/20Sn sample. 



Figure 5: LiOSn sample (eolor online). The eolor eode and dashed box are defined as in Figure]^ 
As for the Lii/20Sn sample the (100) and (010) surfaees are identieal but the strueture has further 
evolved and the Sn atoms are now eompletely expelled from the oxide and reside in ’surfaee’ 
planes on the forming lithia. 
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Figure 6: Fully saturated Li20Sn^ (color online). The same color code as before is used. As the 
lithia saturates the super cell regains an orthogonal shape in the (001) plane, but is skewed along 
the Z-direction. The structure of the (100) and the (010) plane is no longer similar. However, 
the layered structure persists and Sn atoms remain as ’surface’ planes on the lithia. All bonds 
connecting Sn atoms are now longer than the cut off distance, which was defined in Figure]^ 


3.2.3 Li20LSn 

By inserting an additional plane of Li atoms into the LiOSn structure a fully saturated lithia is 
obtained. It is plotted in Figure All the O atoms have now a six-fold coordination and reside at 
octahedral sites where they have five bonds of equal length (1.9 A) to the neighbor Li atoms and 
a single bond to Sn that is 2.3 A. The Li atoms are either two-fold (interface planes) or three-fold 
coordinated (central plane). It is worthwhile noting that none of the Sn-Sn bonds are within the 
defined cutoff distance, as can be seen in Figure]^ The ’surface’ planes therefore appear as weakly 
bound to the enclosed lithia. Along the X-direction the spacing between the Sn atoms is 3.0 A, 
the distance along the Y-direction is 4.9 A, while Sn atoms belonging to two neighbor layers are 
separated by 3.3 A. The shape of the supercell is monoclinic and only the (010) surface is skewed. 
However, the (100) and (010) planes are no longer equivalent. As compared to the pristine SnO 
the volume expansion of the Li20^Sn is 58.59%. This number is consistent with the results for 
Lii/20Sn and LiOSn since they are all composed of a central lithia that is surrounded by planes of 
Sn atoms. 
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4 Discussion 


Our simulation results for the lithiation of SnO align very well with the experimental findings that 
Li behaves as an interealated speeies at low eoneentrations.l^ A similar behavior of Li has been 
observed for Sn02 in a eomputational study by Sensato et As the Li load level inereases four- 
atom clusters form, as shown in Figure [^(b). Due to their strong binding energy these clusters 
retard the emergence of continuous Li stripes. Superstructures made of lines of these low energy 
four-atom clusters should appear instead. 

At higher Li concentrations neighbor Li clusters are more likely to merge and initiate the oxide 
transformation through the proposed ’zipper’ mechanism. Such an evolution is consistent with 
the stripe nucleation experimentally determined for Sn02.^ It cannot be demonstrated that the 
nucleation shown in Figures [^(c) and|^ and the data reported in Reference [^result from the same 
physical / chemical process, but their similarity is striking. However, the oxide transformation is 
likely to be a spontaneous process. 

As the lithia forms and saturates the expelled Sn atoms are equally distributed above and below 
the central O plane. The outcome of this process is that the Sn atoms form ’surface’ planes on 
the layers of the emerging lithia, which can be clearly seen in Figures (c), and Since 
Sn agglomerates into planes outside the lithia, a distribution of various Sn-Sn, Sn-0 and Sn-Li 
bonds occur. This is significantly different from a sample with Sn nano-clusters, which would be 
dominated by Sn-Sn bonds. Our structures with a wide bond distribution could explain the ’exotic’ 
Sn bonds observed in experiments when a SnO sample is lithiated.^^^^^ 

Another important aspect that validates our model is the good qualitative agreement with the 
volume expansion curve measured for LixOSn with respect to SnO and presented in Reference [8l 
As in the experimental data, three distinct regimes can be identified in our computational results 
shown in Figure (i) intercalation at low Li concentrations that results in a slow volume expan¬ 
sion, (ii) oxide transformation when the four atom Li clusters start to merge, causing an abrupt 
volume change, and finally (iii) lithia saturation when the number of Li atoms becomes similar 
to Sn, which slows down the volume expansion. However, the determined Li load levels for the 
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onset of the different regimes do not eoineide with the experimental ones. This diserepancy can be 
explained by the fact that not all the initial SnO layers are lithiated at the same rate. We restricted 
our simulations to a single LixOSn layer and by imposing periodic boundary conditions along the 
X, Y, and Z directions, all the anode layers are assumed to be identical. In reality, as mentioned 
in Reference H parts of the electrode undergo an oxide transformation, while other parts are still 
able to absorb Li as an intercalated species. Hence, the simulated results should be interpreted as 
the upper boundary for the volume expansion during the lithiation process. For example Lii/20Sn 
represents the end point for the oxide transformation in terms of volume expansion and marks the 
transition to the final regime (lithia saturation). The introduction of additional Li further saturates 
the oxide, leading to the formation of LiOSn and finally Li20^Sn. During this phase the slope of 
the volume expansion is reduced, which agrees well with the experimental observations in Refer¬ 
ence m Note that the volume expansion predicted when assuming totally segregated Li20 and Sn 
crystals, the black square in Figure[^ underestimates the experimental value by 7.5 %. In contrast 
the value predicted by our model exceeds the experimental data, but only by 3.1 %. As the dis¬ 
tance between Sn atoms belonging to different layers decreases for increasing Li concentrations 
the volume expansion appears to be caused by the transformation from tin oxide into lithia and the 
formation of Li 20 ^Sn. 

The layered Li20^Sn structure is thermodynamically metastable as compared to a sample 
where Li20 and Sn are segregated. Because the Sn atoms form bonds to Li20^ and the lack of a Sn 
concentration gradient within the planes, the proposed structure should remain sufficiently stable 
as long as the electrode is operated under moderate loading and unloading conditions. The perfor¬ 
mance of a SnO anode has indeed been seen to retain its storage capacity if a reduced voltage range 
is applied.If operated under extreme conditions an agglomeration of Sn into nano-clusters 
might occur for both SnO^^^I^ and Sn 02 .^But this is not the case that has been considered in this 
paper 

A final point supporting the validity of our model is the electron conductivity: in a system 
with completely separated Li20 and Sn the electrons could not efficiently propagate from the 
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metal contact to the active region of the anode or vice-versa, thus limiting the Li reduction or 
oxidation rate. In effect the electron conductivity of oxides is poor and the Sn clusters are not 
necessarily connected to each other to ensure electron transport. Since SnO has been demonstrated 
as a promising anode material, there must be well-conducting electron channels, as for example 
the ’surface’ planes in Figures and which are made of Sn and directly connected to the 
lithia. 


5 Conclusions 

The structure of Lii/20Sn, LiOSn, and Li20Sn has theoretically been determined by lithiating a 
SnO sample. It has been found that during this process the volume expansion undergoes three 
distinct regimes consistent with the experimental observations of Reference [8l (i) slow volume 
changes when Li occupies intercalated sites; (ii) an increased volume expansion when Li reacts 
with O to form lithia; (hi) reduced volume expansion when saturating the lithia. The transformation 
from tin oxide to lithia is likely to be initiated form a LiO stripe resulting from a ’zipper’ nucleation. 

All three phases mentioned above keep the layered and ordered structure of the initial SnO 
crystal. The central LixO layer is embedded between Sn ’surface’ planes. These planes allow us 
to explain the presence of ’exotic’ Sn bonds when SnO is lithiated. The final structure we studied, 
Li 20 ^Sn, is metastable as compared to a fully segregated and lower energy sample with Li 20 
and Sn crystals. In our case the Sn atoms residing between Li20^ appear to be bound to these 
layers, which prevent Sn from segregating and forming nano-clusters if the battery is operated 
under moderate conditions. 

The fact that Sn remains in-between the Li20^ layers challenges the common assumption that 
Sn agglomerates into clusters within a Li20 matrix when SnO is lithiated. This key feature strongly 
influences the structural changes occurring when additional Li atoms are loaded into the SnO 
sample. This issue will be addressed in an upcoming paper. 
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